Modulational instability in a layered Kerr medium: Theory and Experiment 
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We present the first experimental investigation of modulational instability in a layered Kerr 
medium. The particularly interesting and appealing feature of our configuration, consisting of 
alternating glass-air layers, is the piecewise-constant nature of the material properties, which allows 
a theoretical linear stability analysis leading to a Kronig-Penney equation whose forbidden bands 
correspond to the modulationally unstable regimes. We find very good quantitative agreement be- 
tween theoretical, numerical, and experimental diagnostics of the modulational instability. Because 
of the periodicity in the evolution variable arising from the layered medium, there are multiple 
instability regions rather than just one as in the uniform medium. 



PACS numbers: 05.45.Yv, 42.65.Sf, 42.65.Tg, 42.65.-k 

Introduction. The modulational instability (MI) is a 
destabilization mechanism for plane waves. It leads to de- 
localization in momentum space and, in turn, to localiza- 
tion in position space and the formation of solitary- wave 
structures. The MI arises in many physical contexts, in- 
cluding fluid dynamics 0, nonlinear optics 0, plasma 
physics 0, Bose- Einstein condensate (BEC) physics Q, 
and so on. 

The MI was originally analyzed in uniform media, 
mainly in the framework of the nonlinear Schrodinger 
equation. There, MI occurs for a focusing nonlinear- 
ity and long-wavelength perturbations of the pertinent 
plane waves 0, 0, EI 0]- Recently, several experimen- 
tally relevant settings with (temporally and/or spatially) 
nonuniform media have emerged. Such research includes 
the experimental observation of bright matter-wave soli- 
ton trains in BECs induced by the temporal change of 
the interatomic interaction from repulsive to attractive 
through Feshbach resonances [jj. This effective change 
of the nonlinearity from defocusing to focusing leads to 
the onset of MI and the formation of the soliton trains 
@. Soliton trains can also be induced in optical set- 
tings (e.g., in birefringent media Q). Even closer to this 
Letter's theme of periodic nonuniformities is the vast re- 
search on photonic crystals and the experimental ob- 
servations of the MI in spatially periodic optical media 
(waveguide arrays) [jj and BECs confined in optical lat- 
tices |10| . Finally, apart from the aforementioned results 
pertaining to systems that are periodic in the transverse 
dimensions, there exist physically relevant situations for 
which the periodicity is in the evolution variable. Ex- 
amples were initially proposed in the context of optics 
through dispersion management 1 1] , and have since also 
been studied for nonlinearity management both in optics 
[HE3 and BECs 

In this Letter, we present the first experimental real- 
ization of MI in a setting where the nonlinearity is peri- 
odic in the evolution variable, which here is the propaga- 



tion distance. There is a fundamental difference between 
such a periodic setting and a uniform one: In the latter, 
there is a cutoff wavenumber above which MI is not pos- 
sible. In other words, there is a single window of unstable 
wavenumbers. In a periodic medium, however, additional 
instability windows exist for wavenumbers above the first 
cutoff. Our experiments were designed to demonstrate 
this unique feature of the layered structure. In addition 
to our experiments, our investigation includes a linear 
stability analysis |l5j, nj| , which leads to a Hill equation 
(wh ose coefficients are periodic in the evolution variable) 
|17j . The permissible spectral bands of this equation cor- 
respond to modulationally stable wavenumbers and the 
forbidden bands indicate MI. The obtained experimental 
and analytical results are also corroborated by numerical 
simulations. 

Our setup consists of an optical medium with peri- 
odically alternating glass and air layers. The piecewise 
constant nature of the material coefficients leads to a lin- 
ear stability condition (for plane waves) along the lines 
of the Kronig-Penney model of solid state physics 0] 
(generalizations of which with spatially periodic nonlin- 
earity have been considered in |ljj). This allows us to 
compute the MI bands analytically and to compare the 
experimental findings with the theoretical predictions. 

Experimental Setup. In our experiments (see Fig. [T]l. 
an amplified Titanium: Sapphire laser is used to gener- 
ate 150- femtosecond pulses with an energy of 2 mJ at a 
wavelength of A = 800 nm. The beam profile is approx- 
imately Gaussian with a full-width at half-maximum of 
1.5 mm. The laser pulses are split into a pump and a 
reference using a beam splitter (BS1), with most of the 
energy in the pump pulse. After synchronization with a 
variable delay line (DL), the two pulses are recombincd 
at a second beam splitter (BS2) and sent to the peri- 
odic nonlinear medium (NLM) . The reference introduces 
a sinusoidal modulation in the intensity (i.e., an interfer- 
ence pattern), with the period determined by the relative 
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FIG. 1: Experimental setup. BS1 and BS2 are beam splitters, 
DL is a variable delay line, Ml and M2 are mirrors, NLM is 
the layered nonlinear medium, and LI and L2 are lenses. 



angle between the two beams. The angle of the refer- 
ence is carefully tuned by rotating BS2 so that the two 
beams overlap while propagating through the NLM at 
adjustable angles. The NLM consists of six 1 mm thick 
quartz slides separated by air gaps. The glass slides have 
an anti-reflection coating to minimize the loss (the re- 
flection from each interface is 1%). The loss due to back- 
reflections from the slides is included in our numerical 
simulations below, and the effect of double reflections is 
negligible. In our experiments, we used structures with 
air gaps of 2.1 mm and 3.1 mm. The intensity pattern af- 
ter the NLM (at the output face of the last quartz slide) 
is imaged on a CCD camera (Pulnix TM-7EX) using two 
lenses (LI and L2) in a 4-F configuration, with a magni- 
fication of M = 8. The CCD camera captures the central 
region (0.6 mm x 0.8 mm) of the beam. 

The intensity pattern at the output of the NLM is 
recorded both for a high pump intensity (Ip\ = 1.3 x 10 11 
W/cm 2 ) and a low pump intensity (Ip 2 = 9 x 10 8 
W/cm 2 ). In both cases, the intensity of the reference 
beam is 1% of that of the pump. We measure the effect 
of the nonlinearity by comparing the output for high ver- 
sus low intensity. In the latter case, the propagation is 
essentially linear. If the spatial frequency of the modula- 
tion lies inside the instability window, the amplitude of 
the reference wave increases at the expense of the pump. 

Theoretical Setup. Our theoretical model for the beam 
propagation incorporates the dominant dispersive and 
Kerr effects in a nonlinear Schrodinger equation, 



du 
du 



-iv^u _ \ u \ 2u ) < C < 1 (glass) . 



(1) (2) 

ilk V i u -^T)l u | 2u ' ~ 1 <C<L (air), (1) 



2 4 



where space is rescaled by the wavenumber, (£, 77, Q — 
fcW x (x,y, z), and the electric field envelope is rescaled 
using u = (p,2 /n^py-^E. The superscript (j) denotes 
the medium, with j — 1 for glass and j = 2 for air. The 
Kerr coefficients of glass and air are rSp = 3.2 x 10~ 16 



cm 2 /W and 



~ 19 cm 2 /W, respectively. Ad- 
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ditionally, Uq 1 ' = 1.5 and n)^ 1 = 1. The above setting 
(incorporating the transmission losses at each slide) can 
be written compactly as 



du 
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where D(() and N(() are piecewise constant functions in 
consonance with Eq. and 7(C) is the loss rate. Equa- 
tion @ possesses plane wave solutions of the form, 
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where Aq is the initial amplitude. We perform a stability 
analysis by inserting a Fourier-mode decomposition, u = 
uq(() [1 + w(() cos(k^) cos(k v i])] (where w — F + iB is 
a small perturbation) into Eq. J5J. This yields 



d 2 F _ 1 dD dF 
~d~C ~ D~dC~dC 



--k^D 2 + Nk 2 D\u \ 2 



F, (4) 



where k 2 = /c| + k 2 . While one can analyze Eq. Q di- 
rectly, the weak variation of -D(£) can be exploited by 
substituting D(£) with its average and the losses at the 
interfaces can be ignored. [We have checked that this has 
little effect on the results from Eq. . Under these addi- 
tional simplifications, Eq. Q is a Hill equation which for 
the piecewise-constant nonlinearity coefficient is the well- 
known Kronig-Penney model 18] . This can be solved an- 
alytically in both glass and air (with two integration con- 
stants for each type of region) . We match the solutions 
at the glass-air boundaries and obtain matching condi- 
tions at £ = I and C, — L. In so doing, we employ Bloch's 
theorem (and the continuity of F and ^), according to 

which F(CJ = e- ibjC H((), where H is a periodic function 
of period L 18] . This yields a homogeneous 4x4 matrix 
equation, whose solution gives the following equation for 
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+ cos(sif) cos[s 2 (i - I)] = G(k) , (5) 
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uq\ 



and 



k 2 D^(k 2 D^ /4-7V( 2 )|u | 2 ). Therefore, \G(k)\ < 1 im- 
plies stability and \G(k)\ > 1 leads to MI. 

Results. Before discussing our results, it is necessary to 
point out two additional assumptions. First, we assume 
in our numerical simulations that the dynamics is effec- 
tively one-dimensional (ID) along the direction of the 
modulation (i.e., we use k n = and vary k^). Accord- 
ingly, we convert the 2D interference patterns recorded 
on the CCD to ID ones by integrating along the direc- 
tion orthogonal to the modulation. Second, we assume 
that the modulational dynamics of the (weakly decaying) 
central part of the Gaussian beam of the experiment is 
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FIG. 2: (Color online) Comparison of experimental (top), nu- 
merical (middle), and analytical (bottom) results for the 1 
mm glass-2.1 mm air configuration as a function of the di- 
mensionless wavenumber k. For the diagnostics R and \G\ 
(defined in the text), values larger than 1 correspond to MI. 



similar to that of a plane wave with the same intensity. 
We tested both assumptions and confirmed them a pri- 
ori through our experimental and numerical results and 
a posteriori through their quantitative comparison. 

The input field is u — Aq + eoexp(ifc^^), where Aq 
and eo are the amplitudes of the pump and reference 
beams, respectively, and |e | 2 -C |^4o| 2 - For linear prop- 
agation (low intensity, Ip2), the intensity pattern at the 
output of the NLM is approximately the same as that 
at the input; it is about |^4o| 2 + | £ o| 2 + 2A 6q cos(k^). 
For the nonlinear case (high intensity, Ipi), higher har- 
monics are generated and the intensity is about |Ai| 2 + 
2Ai€\ cos(k^) + 2A\C2 cos(2fc^^) + • ■ • , where A% and e n 
{n = 1,2, •••) are, respectively, the amplitudes of the 
pump beam and the nth harmonic at the output of the 
NLM. The Fourier transform (FT) of this latter intensity 
is |Ai| 2 (5(/ 5 ) + Axei5{f t - k ( /2n) + A x ei&{fz + k s /2n) + 
A 1 e 2 6(f i -k ( /ir)+A 1 e25(f ( +k ( /Tr) + - ■ ■ . The ratio of the 
first and zeroth order peaks in the FT is approximately 
equal to the ratio of the amplitudes of the reference and 
pump waves: ro = cq/Aq and r% = e\JA\. (For the ex- 
perimental value of eo = Aq/10, the error introduced by 
this approximation is roughly 1%.) The value of r% in- 
creases with propagation distance as the amplitude of the 
reference increases. In the linear case, ro is constant. We 
use the ratio R — r\ Jtq as a diagnostic measure for both 
our experimental and numerical results (so that R > 1 
indicates growth of the perturbation). This measurement 
is equivalent to the ratio ri(£ = C)/ri(C = 0) (where ( is 
the scaled NLM length) but is more robust experimen- 
tally. In the numerical simulations, the peaks in the FT 
are sharp, whereas they are broader in the experiments. 
Thus, when computing R from the experiments, we used 
the area under the peaks instead of the peak value. 
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FIG. 3: (Color online) Experimental normalized ID intensity 
patterns (top) and Fourier spectra (middle) and numerical 
Fourier spectra (bottom) at the end of propagation (z = 16.5 
mm) of the layered structure with six 1 mm glass slides, 
each pair of which sandwiches 2.1 mm of air. We depict the 
wavenumbers k — 0.0042 (first instability band; left panels) 
and k = 0.017 (second band; right panels). The dashed (solid) 
curves are for high (low) intensity Ipi (Ip2)- 



FigureElshows the ratio R{k) (where k — k^rj^ is the 
sine of the angle between the pump and reference beams) 
for the structure with 2.1 mm air spacings. There are two 
instability bands (quantified by R > 1) within the mea- 
surement range. The appearance of the second band is 
a unique feature of the layered medium that originates 
from the periodicity of the structure in the evolution vari- 
able. The maximum growth of the perturbation in the 
first and second bands appear at k = 6.0 x 10~ 3 and 
k = 1.70 x 10~ 2 , with values of R = 2.05 and R = 1.19, 
respectively. The increase in the modulation is clearly 
visible in the ID intensity patterns (see Fig- El - The po- 
sition of the instability bands is in very good agreement 
with both numerical and theoretical (|G| > 1) predic- 
tions. However, the simulation typically shows a stronger 
instability than the experiment. This results from the 
latter's 3D nature, which is not captured in the simula- 
tion. In the experiment, the spatial and temporal overlap 
of the two beams decreases with increasing fc, resulting 
in weakening of the higher-order peaks. The nonlincar- 
ity also has a lower aggregate strength in the experiment 
because of temporal dispersion. Nevertheless, the simu- 
lations successfully achieve our primary goal of quantita- 
tively capturing the locations of the instability windows. 

The top panels of Fig. show the normalized ID in- 
tensity pattern at the output of the NLM for Ipi (dashed 
curve) and Ipi (solid curve). The left panels are for 
k = 4.2 x 10~ 3 (in the first instability band), and the 
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FIG. 4: (Color online) Same as Fig.[5]but for the 1 mm glass- 
3.1mm air configuration. Here there is a third MI band. 



right ones are for k = 1.70 x 1CP 2 (in the second band). 
The amplitude of the modulation increases for high in- 
tensity cases due to MI. We have also observed in the 
FT of the intensity patterns (middle panels of Fig. 01 the 
appearance of higher spatial harmonics of the initial mod- 
ulation in the instability regions. The first-order peaks 
(the ones closer to k — 0) correspond to the modulation 
of the input beam and are present for both low and high 
intensity. The harmonics correspond to the narrowing of 
the peaks in the spatial interference pattern and appear 
only for high intensity. Strong harmonics are expected 
only within the instability regions. We also observed this 
harmonic generation in numerical simulations (bottom 
panels of Fig. in good agreement with the experi- 
ments. In contrast, we did not observe such harmonics in 
the experiments for k values in the stable regions (R < 1 
in Fig. 2), again in agreement with theory. 

Figure^ shows the instability windows for a structure 
with a different periodicity (with 3.1 mm air between 
each pair of 1 mm glass slides). The instability bands 
shift towards lower k and, as expected by Bloch the- 
ory, the longer spatial period in the structure results 
in a smaller spacing between the instability bands in 
Fourier space. The peaks of the first two bands are at 
k = 5.3 x 10~ 3 and k = 1.46 x 10~ 2 , respectively, and 
a third band appears around k = 2.05 x 10~ 2 . Once 
again, we obtain good agreement between experiment, 
numerics, and theory. It is not straightforward to give 
an intuitive explanation of the precise location/width of 
the stability bands or instability zones [beyond employing 
the simple transcendental expression of Eq. (JSJ)]. How- 
ever, by examining simplified cases (e.g., very narrow air 
gaps between wide glass slides or vice versa) , one can 
extract useful information, such as the decreasing width 
of the MI zones for increasing zone index (which can be 
seen, e.g., in Figs. El and 0. 

Conclusions. We provided the first experimental real- 



ization of modulational instability (MI) in a medium peri- 
odic in the evolution variable. The linear stability analy- 
sis of the pertinent plane waves led to an effective Kronig- 
Penney model, which was used to describe the instability 
bands quantitatively, providing a direct association of the 
MI bands with the latter's forbidden energy zones. One 
of the hallmarks of the periodic medium is the presence of 
additional MI bands as opposed to the single band that 
occurs in uniform media. We found very good agree- 
ment between our theoretical predictions for the modu- 
lationally unstable bands and those obtained experimen- 
tally and numerically. We also observed higher spatial 
harmonics for modulationally unstable beams (another 
characteristic trait of MI). Many interesting extensions 
are possible, including the study of solitary waves that 
result from MI in layered Kerr media. 
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